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Q I Abstract 

In the present paper we numerically construct new charged anti-de Sitter 
black holes coupled to nonlinear Born-Infeld electrodynamics within a certain 
^ ' class of scalar-tensor theories. The properties of the solutions are investigated 

. both numerically and analytically. We also study the thermodynamics of the 

black holes in the canonical ensemble. For large values of the Born-Infeld pa- 
rameter and for a certain interval of the charge values we find the existence of a 
first-order phase transition between small and very large black holes. An unex- 
pected result is that for a certain small charge subinterval two phase transitions 
have been observed, one of zeroth and one of first order. It is important to note 
that such phase transitions are also observed for pure Einstein-Born-Infeld-AdS 
black holes. 
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1 Introduction 



In the last decade anti-de Sitter (AdS) black holes and especially their thermodynamics 
have attracted considerable interest due to the AdS/CFT duality. According to the 
AdS/CFT conjecture the thermodynamics of the AdS black holes is related to the 
thermodynamics of the dual conformal field theory (CFT) residing on the boundary 
of the AdS space [U 12] • In their pioneering work Hawking and Page [5] showed the 
existence of a phase transition between the AdS black hole and thermal AdS space. 
The AdS/CFT duahty provides us with a new tool to study the phase transitions in 
the dual CFT theories on the base of studying the phase transitions of the AdS black 
holes. 

Charged AdS black holes and their thermodynamics have been studied not only 
within the linear Maxwell electrodynamics but also within the framework of the non- 
linear electrodynamics of Born-Infeld. This nonlinear electrodynamics was first intro- 
duced by Born and Infeld in 1934 as an attempt to obtain a finite energy density model 
for the electron [4]. The interest for nonlinear electrodynamics has been later revived 
in the context of string theory. It arises naturally in open strings and D-branes p] 
- [8]. Nonlinear electrodynamics models coupled to gravity have been discussed in 
different aspects (see, for example, [9] - [20] and references therein). In particular, 
black holes coupled to Born-Infeld electrodynamics, including in the presence of a cos- 
mological constant, have been studied in numerous papers [I7],|2T] - |10], and their 
thermodynamics in |?T] - 

The aim of the present paper is to study the charged AdS black holes coupled to 
nonlinear Born-Infeld electrodynamics and their thermodynamics within the frame- 
work of a certain class of scalar-tensor theories (STT). Scalar-tensor theories of gravity 
are among the most natural generalizations of general relativity (GR) and arise in 
string theory and higher dimensional gravity theories [15]. The fundamental question 
is whether objects found and studied in the frame of GR, such as black holes, would 
have different properties in the frame of the scalar-tensor theories. Recent studies 
on asymptotically flat, scalar-tensor black holes coupled to nonlinear electrodynamics 
show that the presence of a scalar field in the gravitational sector leads to interesting 
and serious consequences for the black holes. In the general case the scalar field re- 
stricts the possible causal structures in comparison to the pure Einstein gravity 
For some classes of scalar-tensor theories [IH] there even exist non-unique scalar-tensor 
black hole solutions with the same conserved asymptotic charges. In other words the 
spectrum of asymptotically fiat black hole solutions coupled to nonlinear electrodynam- 
ics in the STT is much richer and much more complicated than in general relativity. 
This shows that the scalar-tensor black holes coupled to nonlinear electrodynamics 
and other sources with non-vanishing trace of the energy-momentum tensor (in four 
dimensions) constitute an area that deserves further study. In view of the current great 
interest in the AdS black holes, the natural next step is to study the scalar-tensor black 
holes coupled to nonlinear electrodynamics in spacetimes with AdS asymptotic. 

In this paper we construct new numerical solutions describing charged AdS scalar- 
tensor black holes coupled to nonlinear Born-Infeld electrodynamics and we make a 
parameter study of their properties. The thermodynamics of the constructed solutions 
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in canonical ensemble is also studied. We found that for large values of the Born-Infeld 
parameter b and for a certain interval of intermediate values of the charge there exist 
a first-order phase transition between small and very large black holes. An interesting 
result is that for a small subinterval one more phase transition exists which is of zeroth 
order since the thermodynamic potential is discontinuous there - it jumps to a lower 
value. Such phase transition of first and zeroth order can be observed also in pure 
Einstein gravity, i.e. when no scalar field is present. 

The phase structure of the Einstein-Born-Infeld black holes and the possible 
Hawking-Page transitions have been examined in [H] and [17]. Hawking- Page phase 
transitions have been studied for different black holes, for example, Einstein-Maxwell 
black holes in spacetime with different dimensions H9] . in three-dimensional space- 
time [501 El], in higher- derivative gravity [521153] . in higher-curvature gravity [M], and 
in unusual electrodynamics which does not restore the Maxwell electrodynamics in the 
weak-field limit [55]. Phase transitions between solitons and black holes in asymptoti- 
cally AdS/Zfc spaces have been studied in [56] . 

2 Formulation of the problem 

The most general form of the action in scalar-tensor theories is 

S = j d'x^ (^Fi^)R - Zi^yg^'d.^d,^ - 2^7($)) + S„,[^rn; 9,.], (1) 

where is the bare gravitational constant and R is the Ricci scalar curvature with 
respect to the spacetime metric g^'^. F{^), Z{^), and U{^) are functions of the 
scalar field $ and their specific choice determines the scalar-tensor theory completely. 
In order for the gravitons to carry positive energy the function must be pos- 

itive > 0), while the non-negativity of the scalar field energy requires that 

2F($)Z($) +3[(iF($)/(i$]2 > 0. The action of the sources is Sm[^rn;gf,u] and there is 
no direct coupling between the sources of gravity and the scalar field in order for the 
weak equivalence principal to be satisfied. 

The action ([T]) is in the so called Jordan frame which is the physical frame but it is 
more convenient to work in the Einstein frame. The relations between the two frames 
are 

g^, = F{^)g^,, (2) 



d^V ( dHFm \' Z($) 
d^J 4 V d^ J 2F{^) ' 

where gf^^ and are correspondingly the metric and the scalar field in the Einstein 
frame. The additional introduction of the functions 

A{^) = (4) 
V{v) = lu{^)F-\^), (5) 
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leads to the following action in the Einstein frame: 



S = -tV / d'^V^ (R - 'Ig^^'d^^du^ - 4V^(<^)) + A\^)g^X (6) 

IottG* j 

where R is the Ricci scalar curvature with respect to the spacetime metric gy,^. In the 
Einstein frame there is direct coupling between the sources of gravity and the scalar 
fields through the coupling function ^(93) and the specific choice of the scalar-tensor 
theory is completely determined by this function and by the potential of the scalar 
field V{ip). 

Here we will consider nonlinear electrodynamics and its action in the Einstein frame 
is given by 



Sm = ^ J d'xy^A\v)L{X, Y), (7) 

where L{X, Y) is the Lagrangian of the nonlinear electrodynamics. The equations 
defining the functions X and Y are 

X = •^^F.^F^'^, (8) 



y = ^^^F,^i^Fr, (9) 

where F^^, is the electromagnetic field strength tensor and ^ stands for the Hodge dual 
with respect to the metric g^^^. 

The Lagrangian of the Born-Infeld nonlinear electrodynamics is 



2b 



X F2 



(10) 



where 6 is a parameter and in the limit 6 — > 00 the linear Maxwell electrodynamics is 
restored. 

The variation of the action ([6]) leads to the following field equations 

= 2dy^d,^ + 2V{v)gy, - 2dxL{X, Y) {p^pFij - ^g^^F^pF"^^ 
-2A\v) [L{X, Y) - YOyHX, Y)] g,,, 



Vy[dxL{X,Y)F^^'' + dYL{X,Y){^Fr''] = 0, (11) 
V.V^if = - Mv)A\v) [L{X, Y) - XdxHX, Y) - YOyHX, Y)] , 

where 

a{^) = (12) 
4 



In the present paper we are interested in spacetimes with anti-de Sitter asymptotic 
structure and thus we will consider the following type of potential: 

V=\k (13) 

where A < is a constant. The coefficient 1/2 is a standard normalizing factor. 

Let us explain the reasons for choosing the potential f[T^ . For this purpose we 
consider the Jordan frame and we require the metric to have AdS asymptotico in this 
frame. Then we have two possibilities - either the factor is finite at infinity 

(i.e. < F($oo) < oo) or it is divergent at infinity. Although the case when -F(<l>) is 
divergent at infinity could be of some interest, it is not generic and in fact is degenerate 
since it corresponds to zero effective gravitational constant at infinity. Here we consider 
only "physically well-behaved" solutions with finit^ F($oo) > . Using the (Jordan 
frame) field equations and especially the equation for $, one can show that the finiteness 
of F($) at infinity requires 

^ v^l d<i> l^F2($) ) ^ 
2Z($)F(<|.)+3(^)'J 

(14) 

lim [/($) = A. 

The simplest (Jordan frame) potential satisfying this conditioijfl, i.e. allowing si- 
multaneously AdS asymptotic for the metric and finite F($oo) > 0, is f/($) = 
AF^($)/F^($oo) which corresponds exactly to the Einstein frame potential ( IT3|) for 
F{^oo) = 1- The conditions (HM written in the Einstein frame take the form 



lim 

>-<I>oc 



2f/($) 



d<S> 



2Z($)F($) + 3(^ 



lim 



lim = 0, lim V{ip) = -A. (15) 

^P^lpOC dip V^VOO 2 

From here it is obvious that V{(p) = A/2 is indeed the simplest choice of a potential 
admitting the desired properties. 

One can show, as we shall see below, that for the scalar-tensor theories under 
consideration, the asymptotic behavior of the factor -F($) is -F($) ~ 1 + const /r^ and 
it guarantees the same anti-de Sitter asymptotic structure simultaneously in both the 
Jordan and the Einstein frames. Moreover, the asymptotic behavior of F{^) guarantees 
that the mass of the Jordan frame solution is the same as that of the Einstein frame 
solution. This can be easily seen from the results in Appendix [Bl 

"'^We mean AdS asymptotic corresponding to a cosmological term A. 
^Throughout this paper, without loss of generahty, we set F{^oo) = 1- 

•^More precisely, we mean the simplest potential in the framework of the class of scalar- tensor 
theories considered in the present paper. For this class of scalar-tensor theories the denominator in 
([14)) is finite. 
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Another reason for considering finite F($oo) comes from the thermodynamics. The 
black hole thermodynamics in scalar-tensor theories is naturally defined in the Einstein 
frame as discussed in Appendix \M When is divergent at infinity the conformal 

transformation does not preserve the AdS asymptotic when we move from the Jordan to 
the Einstein frame. In this case the expected thermodynamic would be rather different 
from that of the AdS black holes if it could be defined at all because of the rather 
unusual spacetime asymptotic. 



3 Basic equations 
3.1 The reduced system 

Here and below we will work in a system of units in which = c = fio/An = 1, where 
Ho is the magnetic permeability of the vacuum. In this system 

[P] = m; [M] = m; [b] = m"^. 

We will consider a static and spherically symmetric spacetime and the ansatz for the 
metric (in the Einstein frame) is 

ds^ = g.^dx^'dx'' = -f(r)e-^^^'Ut^ + ^ + {dO^ + sin^ ed(f)^) . (16) 

f{r) 

An important property of the Born-Infeld electrodynamics is the electric-magnetic 
duality symmetry of the theory [57] - [61]. It is sufficient to study only the mag- 
netically charged case, and the electrically charged solution can be obtained from the 
magnetically charged solution using the electric-magnetic rotation defined by 

{g^,, yp, F^,, P, X, L{X)} ^ {g,,, ^G^,, Q, X, L(X)}, (17) 

where the functions denoted by (..) correspond to the dual solution, P is the magnetic 
charge, Q is the charge of the dual solution, and 



yJiiij — — ^ :z ^ ^ — — 

1^" Q p^lv ' 



dxL{X)^ X. (18) 

In the magnetically charged case the electromagnetic field strength tensor is 

F = Psm9d9 Adc^, (19) 
where P is the magnetic charge. For the functions X and Y we obtain 



^^^^ (20) 



The truncated Born-Infeld Lagrangian is 



L(X)=26(l-\/l + f ). (21) 



Using the metric (fT6!) the field equations reduce to the following system of coupled 
ordinary differential equations 



r - 2fS" - 3f'S' + 2fS'' + - ^fS' 



-4 {V{^) + A\^) [2XdxL{X) - L{X)]} 



(22) 



/" - 2/5" - 3f'6' + 2f6'^ + -/' 



-4 {fv" + V{v) + A\ip) [2XdxL{X) - L{X)]} , (23) 



1 - / - r/' + rf6' = 2r' [V{v) - A\v)L{X)] , 



(24) 



d 



dr 



-dip 



dr 



^KMg-V^ + 4r^e-^a{<f)A\if) [XdxL{X) - L{X)] . (25) 
d(p 



These are four equations for only three unknown functions /, ip and 6 but the self- 
consistency of the system is guaranteed by the Bianchi identity. 

In pure Einstein theory the solution describing Schwarzschild-AdS black holes is 



5E(r) = 0, 
/E(r) = 1 



r 3 



(26) 
(27) 



where Me is the mass of the black hole. The asymptotic of the metric function /(r) at 
infinity in our problem is the same (— |Ar^) and this means that /(r) is unbounded at 
infinity. From a numerical point of view it is convenient to introduce a new unknown 
function m(r) (which is finite at infinity in the considered class of theories) using the 
substitution 



fir) 



2m(r) 



3 



(2^ 



Using the system fl22l) - fl25l) after some manipulations we can obtain the following 
simpler system of ordinary differential equations 



d5 / dip 

dr V dr 



dm 
dr 



dr \ dr J 



-Aa{ip)A\ip) [L - XdxL{X)] - rf 



fdifV 
\dr) 



(29) 
(30) 

(31) 



The first equation fl29|) is decoupled and it can be solved separately once a solution for 
the functions m(r) and ip{r) is obtained. 
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We will consider a class of scalar-tensor theories for which a{(f) = const = a > Qj. 
The investigation of the case with a{(p) = const = a < is similar. Although we 
restrict ourselves to a{ip) = const = a the obtained results are qualitativeljl^ the same 
even for more general coupling functions a{ip) > for which a(0) = a. 

In the chosen ansatz for the metric f|T6l) the temperature is given by the following 
relation 

^ f'{r) e-^« 



4tt 



(32) 

r=rH 



3.2 Asymptotic behavior 

The asymptotic behavior of the functions 6, m and at r — )■ cxd is 

5|™ = ^ + 0(r-«), (33) 
m|,^oo = M - — + -f + 0(r-5), (34) 

¥'|™ = ^ + ^-^^ + 0(r ), (35) 

where C3 is a constant and M is the mass of the black hole in the Einstein frame (for 
the definition of the mass in the AdS spaces see Appendix A). Let us compare the 
asymptotic behavior of these functions in AdS and in asymptotically flat spacetime. 
In AdS spacetime the scalar field decreases as and this is much faster than in 
the asymptotically fiat case where the scalar field decreases as 1/r. The function 6{r) 
also decreases much faster in AdS spacetime, where the leading term is proportional 
to compared to the leading term in the asymptotically fiat spacetime. The 
asymptotic behavior of m(r) up to first order in 1/r is the same. 



3.3 Qualitative properties 

We can obtain some information for the black hole solutions using equations 0251] and 
(12^ and the boundary conditions. We prove that the black holes under consideration 
have simpler causal structure than the black holes in pure Einstein-Born-Infeld the- 
ory in AdS spacetime and we also obtained some information for the behavior of the 
unknown functions. The general properties of the considered scalar-tensor black holes 
coupled to Born-Infeld nonlinear electrodynamics in AdS spacetime can be summarized 
as follows: 

1. The black holes have only one event horizon and extremal black holes do not 
exist. 

2. The scalar field is negative on the event horizon, it has no zeros and increases 
monotonically. 

■^Constant a{(p) corresponds to Brans-Dicke scalar-tensor theory. 

^As the numerical results show the picture is not only qualitatively but also quantitatively close 
for all scalar-tensor theories with the same a(0). 
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3. The metric function S is positive on the event horizon, it has no zeros and de- 
creases monotonically. 

To prove these statements we should first note that it can be easily shown that for the 
Born-Infeld electrodynamics the following inequality is satisfied: 

XdxL{X) - L{X) > 0. (36) 

We will first prove that the considered black holes have only one horizon. Let us 
assume that they have two horizons r_ and r_|_, where r+ > r_. When we integrate 
equation ( l25ll in the interval r G [r_,r+] we obtain 



L{X)] dr > 0. (37) 

The horizons are defined as the points where /(r) = so the left-hand side of the 
equation is zero. But from f p6|) and the fact that we consider a{(p) > it is obvious 
that the right-hand side is strictly positive. So we reached a contradiction and this 
means that two black hole horizons cannot exist. 

In order to prove that an extremal black hole cannot exist we will examine again 
equation ( 125|) . that is 

^ (e-VV^) = ^r'e-'a{^)A\v) [XdxL{X) - L{X)] > 0. 

Let us assume that we have a degenerate event horizon at r^. The right-hand side of 
the equation is strictly positive for every > 0, but the left-hand side is zero and 
we reach a contradiction. 

We can investigate the behavior of the function if{r) in the domain r G [th, oo) 
(jH is the radius of the event horizon) using equation fl25l) . As we already noted the 
right-hand side of (l25l) is positive. Integration of equation ( !25|) from r u io r > th gives 

e-^r^fir)^ = A f r'^e-^ a{ip)A\^) [XdxL{X) - L{X)] dr > 0. (38) 
di^ Jth 

This shows that d^p/dr > in the domain [th, oo). Therefore (/^(r) increases monoton- 
ically. Since the asymptotic value of the scalar field at infinity is zero we can conclude 
that if should be negative on the horizon and monotonically increasing to zero. From 
( 138|) taking the limit r — )■ oo we find that C3 < 0. 

Finally, 5(r) decreases monotonically since the right-hand side of equation fl29|) is 
always negative. The asymptotic value of 5{r) at infinity is zero so we can conclude 
that it is positive on the event horizon and monotonically decreasing to zero. 

All of the above presented propositions remain valid for different nonlinear electro- 
dynamics for which the relation (!36l) holds. 

^We are not interested in the case = because this corresponds to naked singularity. 



dr 



dr 



I + 

4 j r^e-^a{^)A\^) [XdxL{X) - 
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3.4 Dimensionless quantities 

The reduced field equations (!29l) - (!3Tl) are invariant under the rigid rescahng 

r ^ Ar, m ^ Am, A A^^A, b \-%, P XP. (39) 

Therefore, one may generate in this way a family of new solutions. Thus, the mass, 
the temperature and the entropy of the new solutions are given by the formulas 

M^AM, T^A-^T, S^X'^S. (40) 

As a consequence of the rigid symmetry we may restrict our study to the case 
A = — 1 by introducing the dimensionless quantities 



r I, 1,72 n ^ 

r y, m y, bl , P y. 



M ^ ^, T^IT, S^^, (41) 



where the choice of A = — 1 requires 



3.5 Posing the Boundary- Value Problem 

The domain of integration is half-infinite r G [th, oo) where the horizon, which is the 
left boundary, is a priori unknown. We impose the following boundary conditions: 

— At infinity: 

lim m(r) = M, (43) 

r— )-cxD 

lim 6{r) = lim (p{r) = 0, (44) 

r— >-oo r— >-oo 

where M is the mass of the black hole in the Einstein frame (see Appendix A); 

— On the horizon: 

firn) = 0. (45) 
The following regularization condition should also be fulfilled on the horizon: 



{M^)A\^)[XdxL{X) - L{X)]} 



0. (46) 



r=rH 



In this way the governing equations ( 129|) - (13T]) together with the above boundary 
conditions (143|) - (H6|) . complement a well-posed boundary value problem. Since the 
position of the event horizon (the left boundary of the integration interval) is not fixed, 
moreover it is unknowrill and there is no analytic way to determine it, we can use one of 



^Such BVPs are known in mathematical physics as BVPs of the Stefan kind. 
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these boundary conditions to determine it. We can introduce a new shifted independent 
variable x = r — th- In this way the domain of integration is completely defined - 
X G [0, oo) and the radius of the horizon will appear explicitly in the equations as an 
unknown parameter. So, we have one parametric problem, which can be considered as 
a spectral-like nonlinear BVP. 

Let us emphasize that there is an alternative formulation of the problem where the 
horizon is an input parameter. In this case we do not need a shifting map and the mass 
M becomes the sought parameter. We used this method of solution for justification of 
our results. 

3.6 Method of solution 

The system of ordinary differential equations fl29l) - fl31l) is coupled and nonlinear. It 
does not admit a global analytic solution and only some local expansions are possible. 
This is the reason to solve it numerically by using the reliable Continuous Analog of 
Newton's Method, which is distinguished by its quadratic convergence in the vicinity 
of a localized root [62l l63][^ . Yet, the method has been already successfully applied to 
similar problems (see, for example, [HU [381 EHl US]). 

4 Discussion of numerical results 

As we noted in the previous sections we will consider constant parameter a > (a is 
defined by equation ( IT2i) ). The results do not change qualitatively for different a > 
and we will present the results for a = 0.01. 

In Figs. [1] and [2] the metric functions m(r) and 6{r) and the scalar field ip{r) for 
P = 0.07, b = 10, a = 0.01 and for several values of the black hole mass are presented. 
An important justification of our results is that, up to the leading terms, the asymptotic 
form of the numerically obtained functions S{r), m{r), and ^{r) at infinity is the same 
as what is predicted in fl33|) . flMl) and fl35|) . respectively. 

In Fig. [3] the horizon and the temperature T as functions of the mass M for 
sequences of black hole solutions are plotted for 6 = 10 and a = 0.01. Fig. H] represents 
two magnifications of the enclosed regions on the T{M) plot in Fig. |3l In Fig. [5] the 
value of the scalar field on the event horizon and the constant which occurs 
in the asymptotic expansion of the function ip{r) at infinity (see equation (15^ ) as 
functions of the mass are plotted. In Figs. [6] and [7] the quantities th, T, ipH and C3 as 
functions of the charge P are plotted. 

When hP^ > 1/8 extremal black hole solutions exist in Einstein-Born-Infeld theory 
in AdS spacetime. In the presence of scalar field extremal black holes cannot exist as 
we have already proved. As we see from the T(M) plot for P = 0.113 and 6 = 10 (this 
combination satisfies bP^ > 1/8) in the right panel of Fig. [3l when M decreases the 
temperature first reaches values close to zero and afterward increases. The reason is 
that the absolute value of the scalar field rises considerably in this region, as we can 

^Let us note that in our case the roots appear to be functions in Banach space. 
^This value is consistent with the current observational data |45) . 
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see in Fig. |5l and prevents the reaching of an extremal black hole. The same situation 
can be seen from the plot T{P) in Fig. E]when we increase the magnetic charge P (see, 
for example, the curve corresponding to M = 0.10). 

The numerical results are obviously subject to and possess the properties in Sub- 
section 13.31 All of the dependences have similar qualitative behavior in the parameter 
space that we have studied. The only exception is the function T(M) because its qual- 
itative behavior depends also on the value of the parameter b (more details on this 
subject will be given in the next sections). 
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Figure 1: The metric functions 6{r) (left panel) and m(r) (right panel) for several 
values of the black hole mass and for P = 0.07, b = 10, and a = 0.01. 



0.000 




Figure 2: The scalar field ip{r) for the same values of the parameters as in Fig. [TJ 
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Figure 3: The radius of the black hole horizon (left panel) and the temperature T 
(right panel) as functions of the mass M for sequences of black hole solutions (6 = 10 
and a = 0.01). Note that on the right panel we denoted "Enclosed region 1" and 
"Enclosed region 2" which are expanded in Fig. HI 




Figure 4: Magnifications of "Enclosed region 1" (left panel) and "Enclosed region 2" 
(right panel) of the T(M) plot in Fig. H 
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0.30 



Figure 5: The scalar field on the horizon ipn (left panel) and the constant C3 (right 
panel) in the asymptotic expansion of the function ip{r) at infinity (see equation (l35l)). 
as functions of the mass M for sequences of black hole solutions (6 = 10 and a = 0.01). 




Figure 6: The radius of the black hole horizon th (left panel) and the temperature T 
(right panel) as functions of the charge P for sequences of black hole solutions {b = 10 
and a = 0.01). 
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Figure 7: The scalar field on the horizon (pH (left panel) and the constant C3 (right 
panel) in the asymptotic expansion of the function ip{r) at infinity (see equation (!35|l ). 
as functions of the charge P for sequences of black hole solutions (6 = 10 and a = 0.01). 
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5 Thermodynamic phase structure 

To study the thermodynamic^ of the black holes under consideration we need to find 
the conserved charges of our system first. Throughout the years many methods have 
been developed to calculate the conserved charges of the gravitational configurations. 
A related problem is the computation of the gravitational action of a noncompact 
spacetime. When evaluated on noncompact solutions, the bulk Einstein term and 
the boundary Gibbons- Hawking term are both divergent. The remedy is to consider 
these quantities relative to those associated with some background reference spacetime, 
whose boundary at infinity has the same induced metric as that of the original space- 
time. This substraction procedure is however connected to difficulties: the choice of 
the reference background is by no means unique and it is not always possible to embed 
a boundary with a given induced metric into the reference background. A new method, 
free of the aforementioned difficulties, was proposed in [65]. This method (called the 
counterterm method) consists of adding a (counter)term to the boundary at infinity, 
which is a functional only of the curvature invariants of the induced metric on the 
boundary. Unlike the substraction procedure, this method is intrinsic to the spacetime 
of interest and is unambiguous once the counterterm which cancels the divergences 
is specified. In the present work we use, namely, the counterterm method in order 
to compute the mass and gravitational action. The technical details are presented in 
Appendix [Bl 

We will study the thermodynamic phase structure for black holes in the canoni- 
cal ensemble where the magnetic charge P is fixed. The associated thermodynamical 
potential in the case of the canonical ensemble, as one should expect from pure ther- 
modynamical considerations, is the Helmholtz free energy -F(T, P) = M — TS. This 
can be proven with more rigor by calculating the Euclidean action as it is done in 
Appendix [B] . 

Let us note that our definition of the free energy is intrinsic. Our computation 
based on the counterterm method makes no reference to any other solution of the 
field equations. This should be contrasted with the method used in ^7] where the 
authors compute the action using the extremal solution (with fixed charge) as reference 
background. As we already showed in our case there are no extremal solutions and 
therefore the reference background method of 0^ is unapplicable. 

The thermodynamic phase structure depends on the value of magnetic charge P 
and the parameter b in the Born-Infeld Lagrangian. The numerical results show that 
for large values of the parameter b the function T = T{rH, P) has two infiection points 
defined by 

These infiection points, which we will denote by P^}^ and P^}^-, separate three charge 
intervals given by P < P^ll, P^ll < P < Pj2 and < P- Below we will consider 
each of these intervals in detail. 



^°Some clarifying comments on the black hole thermodynamics in scalar-tensor theories are given 
in Appendix VK\ 
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For small values of the parameter b the inflection points disappear and the phase 
structure is the same for all the charges. The numerical results show that for all of the 
studied values of the parameter a (including the case a = which corresponds to pure 
Einstein gravity) the critical value of the parameter b which separates the regions of 
small and large b, is focrit ~ 0.6. 

5.1 Phase structure for large values of the parameter b > 6crit 

As we noted for large values of the parameter b > 6crit the function T = T{rH, P) 
has two inflection points - P^^^^j. and P^'^l, which gives us three charge intervals defined 

by P < Peril) Pciit < P < Peril Peril < P- Formally we call these three cases 
black hole with small, intermediate and large charge respectively. On its own side the 
intermediate charge interval has to be divided into three subintervals Peril < P < Ppl^ , 

PjSi^ < P < P^^ and Pp^h^ < P < P^ril (Pph and P^^ are defined and discussed in 
details in the following subsections). As we show below, there are phase transitions for 
charges belonging to the second and the third subintervals. 

All the results for the scalar-tensor black holes given in this subsection are computed 
for fixed values of the parameters b = 10 and a = 0.01. For this choice of parameters 
the values of the critical charges are 

= 0.09397 and Pe^.^jl = 0.2901. The value of 
Pej.jt is the same (within the numerical error) as in the pure Einstein theory (i.e. for 
a = 0). The reason is that the corresponding inflection point, defined by equation (H71) . 
is situated in a region of the parameter space where the contribution of the scalar field 
is very small. The value of Peril in the case of pure Einstein theory - P^^^^l = 0.09395, 
is a bit smaller than in the case when scalar field is present. The numerical results for 
larger a confirm that Peril increases slightly when we increase the value of the coupling 
constant a. 

The values of the charges Pp^^ and Pp^^ for the chosen values of the parameters 
b = 10 and a = 0.01, are Pp^-* = 0.1087 and Pp^'' = 0.1093 and the corresponding values 
in pure Einstein theory are Pph^„=o ~ 

0.1085 and P^l\^o = 0.1091. The numerical 

experiments show that the values of Pp^'' and Pp^^ increase slightly with the increase 
of a. 

The temperature T as a function of the horizon radius th is shown in Fig. [H] for the 
cases of small (for example P=0.07), intermediate (for example P=0.115), and large 
(for example P=0.40) charges P. The free energy P as a function of the temperature 
T for small and intermediate charges is shown in Fig. [9] and magnifications of specific 
regions are presented in Fig. [TOl The function P(T) for various charges, including large 
charges, is given in the left panel of Fig. [T21 In Figs. [TT] and [12] (right panel) we also 
give the dependences T{rH) and P(T) for Born-Infeld-AdS (BIAdS) black holes in the 
pure Einstein gravity for comparison. The figures show that the thermodynamic phase 
structures for the scalar-tensor BIAdS black holes and the Einstein BIAdS black holes 
are qualitatively the same except in the case of charges satisfying bP^ > 1/8 where 
additional phases of the scalar-tensor BIAdS black holes are present. Below we consider 
in detail the thermodynamic phase structure for the intervals of small, intermediate 
and large charges. 
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Figure 8: The temperature T as a function of the radius of the horizon rn for sequences 
of black hole solutions {b = 10 and a = 0.01). 



The phase structure and the phase transitions can be alternatively studied with 
the so-called off-shell formalism or, in other words, off-equilibrium formalism. It is 
thoroughly presented in ^66j and applied, for example, in the papers on black-hole phase 
transitions by Myung and coauthors [501 ED HZ]- The advantage of that formalism 
is that it gives a simple and clear interpretation of the thermodynamic phases. In 
the canonical ensemble the different phases of the system, both stable and unstable, 
correspond to extrema of the so-called off-shell free energy. In that formalism, the origin 
of the branches in the T — F diagram, that correspond to the different phases, can be 
also simply explained. The phase transitions have been studied with the application 
of the off-shell formalism in Appendix O 



5.1.1 Phase structure for small charges 

In the region of small charges P < -Pcrit, as one can see from Fig. [131 fo^' arbitrary 
temperature T > Tq there are two branches of solutions for the black hole radius. 
These two branches correspond to small and large black holes, respectively. The small 
black holes have negative specific heat Cp which is defined according to the formula 



Cp = T ^ 
^ dT 



dS „ „ f dT^'^ 



2TxTrH 



\drh 



(4^ 



and therefore they are unstable. The large black holes have positive specific heat and 
consequently they are locally stable. Moreover, as can be seen in Fig. [T3] (right panel) 
the free energy of the large black holes is smaller than the free energy of the small 
black holes and therefore the large black holes are thermodynamically favorable and 
dominate the thermodynamic ensemble. 
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Figure 9: The free energy F as a function of the temperature T for sequences of black 
hole solutions {h = 10 and a = 0.01). Note that in the figure we defined "Enclosed 
region 1" and "Enclosed region 2" which are expanded in Fig. nH] 




Figure 10: Magnifications of the enclosed regions of the F{T) plot in Fig. [91 The left 
panel corresponds to "Enclosed region 1" and the right panel corresponds to "Enclosed 
region 2". 
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Figure 11: The functions T^th) and F{T) for sequences of black hole solutions in pure 
Einstein theory (6 = 10 and a = 0.0). 
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/M).300 
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Figure 12: The function F{T) for sequences of black hole solutions with various charges 
P, including large charges, in the case when a — 0.01 (left panel) and for pure Einstein 
theory - a = 0.0 (right panel). 



19 




— ] 1 1 1 1 1 1 1 1 1 1 1 -0.05 -I h 1 1 1 1 1 1 1 1 

0.0 0.4 0.8 1.2 1.6 2.0 0.2 0.4 0.6 0.8 1.0 



Figure 13: The functions T(ri^) and F{T) for small charges where 6 = 10 and a = 0.01. 
The curve corresponding to the critical charge pI:^^ is denoted by the dotted line. 



5.1.2 Phase structure for intermediate charges 

In the intermediate charge interval defined by Pj^^j^ < P < P^'^l, there are four branches, 
as shown in Figs. [HHTTl These branches will be called branches of very small (VSBH), 
small (SBH), large (LBH), and very large (VLBH) black holes, respectively. 

Below we describe the phase structure in the subintervals pI:^1 < P < P^)^ , 

(1) (2) (2) (2) 

Ppj^ < P < Ppj^ and Ppj^ < P < Pcrit- The differences between the three subintervals 
are that in the first subinterval no phase transitions exist, in the second subinterval 
phase transitions of first and zeroth order exist, and in the third subinterval only 
first order phase transition exists. Such subintervals with the corresponding phase 
transitions are also present for pure Einstein-Born-Infeld-AdS black holes. 



Intermediate charge subinterval P^^.-^ < P < P^^ 

In this subinterval Tq < T^-m 0- The Tirn) and F{T) dependences are shown in 
Figs, [m and US] and the phase structure is described as follows. For low temperature 
T satisfying Tq < T < Tmin there are two solutions for the black hole horizon radius 
corresponding to LBH and VLBH. From the dependence T(r//) shown in Fig. [T^one 
can conclude that the LBH are thermodynamically unstable since they have negative 
specific heat, while the VLBH have positive specific heat and therefore they are locally 
stable. At T = Tmin the origin of two new branches of solutions appears and for 
^min < T < T^Q^y. we have four black holes with the same temperature and different 
radii - VSBH, SBH, LBH, and VLBH. The VSBH and LBH are unstable since their 
specific heat is negative, while the SBH and VLBH have positive specific heat and 
consequently they are locally stable. At T = T^ax two of the branches coalesce and 
disappear, leaving only two branches for T > T^ax corresponding to unstable VSBH 

^"'^By Tq and Ty^in we denote the two minima of the T{rH) function as can be seen in Fig. [M] 
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Figure 14: The function T(r//) for intermediate charge pI:^^ < P < P^^i (P = 0.107, 
6 = 10 and a = 0.01). 
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Figure 15: The function F{T) for intermediate charge p£\| < P < P^^l^ {P = 0.107, 



6 = 10 and a = 0.01). The right panel is a magnification of the enclosed region of the 
left figure. 



and locally stable VLBH. From T = Tq to arbitrary large temperature the free energy of 
the VLBH, corresponding to the leftmost branch in Fig. [15] (left panel) (corresponding 
to the rightmost branch in Fig. [T^ . is the smallest which means that for arbitrary 
temperature the VLBH dominate the thermodynamic ensemble. 

Intermediate charge subinterval P^l^ < P < P^f^ 

In this subinterval To < T^m- The T^rn) and F(T) dependences are shown in 
Fig. [16] and the phase structure is described as follows. For Tq < T < Tmin we have 
two branches - unstable LBH and locally stable VLBH. At T = Tmin two new branches 
appear and for Tmin < T < Tmax there are four branches corresponding to unstable 
VSBH, locally stable SBH, unstable LBH and locally stable VLBH. At Tmax two of the 
branches coalesce and disappear and we are left with only two branches of unstable 
VSBH and locally stable VLBH. From Tq to Tmin the VLBH have the smallest free 
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energy and they dominate the thermodynamical ensemble. At T = Tmin the situation 
changes and for temperatures in the interval from T^i^ to Tph the ensemble is dominated 
by the SBH which have the smallest free energy. At this point T = Tmin phase transition 
between VLBH and SBH occurs and it is zeroth order phase transition since the free 
energy is discontinuous at that point - the free energy jumps to a lower valu^. Such 
unusual phase transitions of zeroth order have been observed for example in [67 | [68| [69] . 

As we said for temperatures in the interval Tmin < T < Tph the ensemble is dom- 
inated by the SBH. As we can see in Fig. [16] for arbitrary T > Tph the ensemble is 
again dominated by the VLBH which have the smallest free energy. This means that 
at the point T = Tph we observe a phase transition between the SBH and VLBH and 
it is classified as a first order phase transition since the free energy is continuous at 
T = Tph and the first derivative of the free energy is discontinuous at that point. 

Intermediate charge subinterval Pj^^ < P < Pj^^^l 

In this subinterval Tmm < Tq. The T{rH) and F{T) dependences are shown in 
Fig. [T7| and the phase structure is described as follows. For Tmin < T < Tq there are 
two branches corresponding to unstable VSBH and locally stable SBH. At T = Tq 
two new branches appear and for temperatures Tq < T < T^ax we have four branches 
corresponding to unstable VSBH, locally stable SBH, unstable LBH and locally stable 
VLBH. At T = Tjnax the branches of SBH and LBH disappear and we are left with 
unstable VSBH and locally stable VLBH for T > Tmax- From Tmin to temperature Tph 
the free energy of the SBH is the smallest and the SBH dominate the thermodynamic 
ensemble. At temperature T = Tph the free energy of the VLBH becomes smaller than 
the free energy of the SBH and for all temperatures T > Tph the VLBH dominate the 
thermodynamic ensemble. At the point T = Tph we have a first-order phase transition 
from SBH to VLBH. 

""^^We have studied the zeroth-order phase transition observed at that point with the off-shell for- 
malism in Appendix [C] 
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Figure 17: The functions Tiju) and F{T) for intermediate charge P^^ < P < 
(P = 0.113, 6 = 10 and a = 0.01). 




Figure 18: The functions T{rH) and F{T) for large charge P^i < P {P = 0.35, b = 10 
and a = 0.01). 



5.1.3 Phase structure for large charges 

In the case of large charges defined by P^j.^^ < P, there are two branches - small and 
large black holes, as Fig. [18] shows. The small black holes have negative specific heat 
and consequently they are unstable. The large black holes are locally stable since their 
specific heat is positive. For arbitrary temperature the free energy of the large black 
holes is smaller than the free energy of the small black holes. Therefore the large black 
holes dominate the thermodynamic ensemble. 



5.2 Phase structure for small values of the parameter b < bent 

The phase structure for small values of the parameter b is much simpler because the 
function T = T{rH, P) has no inflection point. Therefore the phase structure is qualita- 
tively the same for all values of the magnetic charge P. The functions T{rH) and F{T) 
are plotted in Fig. [T9]for b = 0.1 and a = 0.01. Only two branches exist corresponding 
to unstable small black holes and locally stable large black hole. The large black holes 
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Figure 19: The functions T{rH) and F{T) for a small value of the Born-Infeld parameter 
b = 0.1 (a = 0.01). 




Figure 20: The functions Tirn) and F{T) for sequences of black hole solutions in pure 
Einstein gravity (a = 0.0) for a small value of the Born-Infeld parameter h = 0.1. 



have lower free energy and therefore they dominate the thermodynamic ensemble. 

The functions T{rH) and F{T) for 6 = 0.1 are shown in Fig. [20] in the case of pure 
Einstein theory which corresponds to a = 0.0. This case is qualitatively the same as 
the case when the scalar field is nontrivial except for large charges which fulfill the 
inequality bP"^ > 1/8. For these values of b and P, when a = 0.0, the branch of 
small black holes disappear because of the occurrence of an extremal black hole (see 
for example the curve corresponding to P = 1.125 in Fig. |20|) . 



6 Conclusion 

In the current paper we have constructed numerically new scalar-tensor black hole 
solutions in AdS space-time coupled to Born-Infeld nonlinear electrodynamics. The 
properties of the solutions have been thoroughly studied through a combination of 
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analytical and numerical techniques. The black hole solutions in the considered class 
of scalar-tensor theories are determined uniquely by the asymptotic charges of the 
black holes - in our case the mass and the magnetic charge. An important property is 
that in comparison to the corresponding solutions in General Relativity, namely, the 
Einstein-Born-Infeld-AdS black holes, the black holes considered here have a simpler 
causal structure - they have neither inner nor degenerate horizons. 

The thermodynamics and the phase structure in the canonical ensemble of the 
obtained black holes have also been studied. Their phase structure is similar to that 
of the Einstein-Born-Infeld-AdS black holes in GR and the difference is manifested at 
large charges. This means that in the case under consideration there are no extremal 
black hole solutions which actually leads to one additional branch of solutions in the 
thermodynamical phase-space structure for large values of the charge. 

The thermodynamical phase-space structure depends on the Born-Infeld parameter 
b and the magnetic charge P. A critical value 6crit ~ 0.6 has been found numerically. 
For b < 6crit the phase diagram is the same for all values of the magnetic charge P. 
There are two branches corresponding to unstable, small black holes and locally stable, 
large black holes. In this case no phase transitions are possible. 

In the case of 6 > 6crit) the structure of the phase diagram changes with the variation 
of the magnetic charge P and three charge intervals are observed - small, intermediate, 
and large charges. For small and large values of the magnetic charge no phase transi- 
tions are possible and the phase diagram has only two branches which corresponds to 
two phases: unstable, small black holes and locally stable, large black holes. The in- 
termediate charges are the most interesting because phase transitions can exist in this 
case. The phases are four: VSBH, SBH, LBH, and VLBH. Among them two are locally 
stable - the SBH and the VLBH - and phase transitions between the two black holes 

(2) (2) 

are possible in some cases. In the intermediate charge subinterval P^^ < P < P(,j.j(. 
first order phase transition between SBH and VLBH exists - for small values of the 
temperature (Tq < T < Tph) the SBH are thermodynamically favorable because they 
have lower free energy. When we increase the temperature we reach a certain value of 
the temperature Tph where phase transition to VLBH occurs and for all temperatures 
above that value T > Tph the VLBH arc thermodynamically favorable. This phase 
transition is of first order because the first derivative of the free energy is discontinu- 
ous at that point. An important observation is that for a small intermediate charge 
subinterval P^l^ < P < P^^^ two phase transitions exist. One of them is the first order 
phase transition we just described and the other one is the zeroth order phase transi- 
tion. In this case for small values of the temperature (Tq < T < Tmin) the only locally 
stable branch is the VLBH. At temperature Tj^in another locally stable branch appears 
- SBH, which have the smallest free energy and dominate the thermodynamical en- 
semble. This means that for T — Tmin phase transition between VLBH and SBH exists 
and it is of zeroth order since the free energy is discontinuous at that point. When 
we increase the temperature the point T = Tph is reached where the first order phase 
transition between SBH and VLBH, that was described above, occurs. For T > Tph 
again the VLBH dominate the thermodynamical ensemble. For all the other interme- 
diate charges phase transitions are not possible since the free energy of the VLBH is 
always the smallest. Similar zeroth and first order phase transitions for certain charge 
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intervals are present also in the corresponding solutions in General Relativity. 

Let us finish with prospects for future work. As in the asymptotically fiat case 
[ISj, in AdS spacetimes the scalar-tensor theories with a{ip) = (3ip {/3 < 0) coupled to 
nonlinear electrodynamics admit non-unique black hole solutions [70]. From AdS/CFT 
perspectives these non-unique black hole solutions could be considered as second order 
phase transitions in the dual field theory living on the boundary. Detailed considera- 
tions will be given in a forthcoming paper [70]. Phase transitions due to non-unique 
AdS black hole solutions in Einstein-Maxwell-dilaton gravity with appropriate dilaton 
coupling function have been obtained quite recently in [71] • Let us note however that 
the situation in the case of scalar-tensor AdS black holes coupled to the nonlinear 
electrodynamics is rather different in comparison with the Einstein-Maxwell-dilaton 
gravity. 
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A Black hole thermodynamics in scalar-tensor 
theories 

In this appendix we would like to make some comments on the thermodynamics of the 
black holes in scalar-tensor theories. The definitions of entropy, temperature, and mass 
(energy) of the scalar-tensor black hole are naturally related to the Einstein frame. The 
reasons are the following. The Einstein-frame area of the horizon can be interpreted 
as entropy of the horizon since it is non- decreasing in classical physical processes. The 
entropy of the event horizon in the Einstein frame Se, for example, is proportional to 
one fourth of the area of the horizon. In the Jordan frame, however, the area of the 
event horizon may decrease in classical processes [72] and a generalized definition of 
entropy Sj needs to be introduced. The proper definition is [721 1121 [HI 175] 



This definition is also consistent with the Euclidean and the Noether charge method 
Using relation ([2]) we find that 




(49) 




(50) 
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the definitions of the entropy in both conformal frames are in agreement. In the last 
two equations quantities ^^^^ and ^"^^g are the determinants of the induced metrics on 
the horizon in the Jordan and in the Einstein frames, respectively. 

Further, it can be proved that the standard definitions of black-hole temperature 
are invariant with respect to conformal transformations that are regular on the event 
horizon [76]. Since the conformal transformations (j2]) in the present paper are regular 
on the horizon (and everywhere outside the horizon), the temperature in the Einstein 
frame coincides with the temperature in the Jordan frame. 

The measure of the internal energy of compact objects in STT needs also to be 
specified. In the general case the masses of the black holes in the Einstein frame and in 
the Jordan frame do not coincide. The mass in the Einstein frame is positive definite, 
can only decrease by emission of gravitational (and scalar waves), and possesses other 
natural energy-like properties unlike the mass in the Jordan frame. In other words 
the mass in the Einstein frame is the true energy of the compact objects within the 
framework of the STT. For more details on the mass of compact objects in the STT 
we refer the reader to [771 ES EH EOl [381 ES]- For the particular solutions studied in 
the current paper, however, the scalar field decreases fast enough and the black-hole 
masses in both conformal frames are the same. 



B Derivation of the free energy 

The action with surface terms and counterterms taken into account is given by 
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(51) 



where 7^,^ is the metric induced on the boundary dAi and is the trace of the extrinsic 
curvature of dAi while Sctll] is the counterterm action. We have also defined £ = 
a/— 3/A and the boundary stress-energy tensor is given by 
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(52) 



If the boundary geometry has an isometry generated by the Killing vector K, T^^K^ 
is divergence free which gives the conserved quantity 



Qk = / da>'%,K'' (53) 

associated with the closed surface a in the boundary dAi. In the case when K = d/dt, 
the conserved quantity is the mass. 

For spacetimes with AdS asymptotic the following simple counterterm was proposed 
in [65] (see also [8T] ) 
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(54) 



where TVy^] is the Ricci scalar curvature of the boundary metric 7^,^. With this coun- 
terterm for the asymptotic metric 

ds'^-(l- — + '^) de + '^l^ + {de' + sin^ ed<p') (55) 
we have up to the leading order 



SttTu = ^ + 0{r-'') (56) 



which gives the mass 



= M. (57) 

In other words this shows that the limit lim^^-oo ''^('") = M of the metric func- 
tion m(r) (see eq. (l28|l ) has to be identified with the mass of the solutions with AdS 
asymptotic. 

Our next task is to calculate the Euclidean action. For definiteness we shall consider 
the electrically charged case. First to make the action Euclidean, the time coordinate 
should be made imaginary by substitute t = ir. This makes the metric positively 
definite, namely 



ds'' = f{r)e-^'^'-^dr^ + ^ + [dO^ + sin^ Odcj)'') . (58) 

In order to eliminate the conical singularity at the horizon r = rn, the Euclidean 
time coordinate r has to be periodic with period /3 = 1/T where T is the Hawking 
temperature associated with the black hole horizon. The very Euclidean action is given 
by 



/ = /bulk + I0 + let + Id (59) 

where 



'bulk 



(60) 
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/e = / d'x^Qi^l (61) 
/cf = ^ / f/=^xy^n^A,F'^'^^^((^)axL. (63) 

471" JdM 

The bulk integral is on compact region Ai with boundary dAi . The boundary dAi 
lays at finite radius r = vb and has topology x S*^ where 5"^ stands for the periodic 
Euclidean time coordinate r. The last term Id is a charge fixing term, in other words 
we have added this term since we consider a fixed charge ensemble. 

In order to compute the Euclidean action and to show that I/f3 = M — TS we will 
need an auxiliary tool, namely the Komar integral 

Mb = -^ [ V^CdS^, (64) 



^vr J SI 



where 5*^ is a sphere with radius rs and ^ is the timelike Killing vector. Obviously the 
Komar integral is divergent for — > oo but it is finite for finite r^. More precisely, 
for the AdS asymptotic (l55ll we have 

Mb = M+^ + 0(— ). (65) 

Using Stokes theorem we can present the Komar integral as a sum of integral on 
the horizon % and a bulk integral 

Mb = -^ [ V^CdS,, [ R,,Cd^^ (66) 

where S is the 3-dimensional space bounded by S% and "H. 
It is well known that the integral on the horizon gives just 

_ _L /■ V^'^-dS^, = 2TS (67) 



H 



where T is the Hawking temperature and S is the black hole entropy which is one 
fourth of the horizon area. On the other hand using the field equations we find 

R,uC = - 2dxLF^aFJ'C + 2^'(<^) [2XdxL - L] (68) 
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Now we introduce the electric potential $ by the equation = —^"^F^a and taking 
into account the field equations for the electromagnetic field, the second term on the 
right-hand side can be presented in the form 



- 2dxLF^aF,''C = 2dxLF^d^^ = 2V„ {^dxLF^) . (69) 
In this way we obtain 

+ ^ ! A\^) [2XdxL-L]^^dW. 



(70) 



With the help of the Stokes theorem we find 



--^ / V, {<^dxLFn dT.^ = I ^^dxLF^'^dS,^ + / <!>dxLF^^''dS^^ = 
0{—) + / dxLF^^'^dS,^ = 0{—) + 2^nQ (71) 

where we have taken into account that at large we have ^{tb) ~ and that the 
electric potential is constant on the horizon. Also, the electric charge of the black hole 
is given by 

Q = ^ [ dxLF'^'^dS^^. (72) 



Summarizing, we obtained the following relatioij^ 



Mb = 2TS + 2$^Q ' ^ ^/^^^^ '^'^^^ ^'^^^^^ ' ^] ^^^^^ (73) 

which holds up to terms of order 0{r^\ Having this important relation we can go 
back to computation of the Euclidean action. Using the field equations the bulk term 
can be presented in the form 



/buik = -^/ Vgd'x-^f A\ip){2XdxL-L)^d^x. (74) 
8vr Jm 47r 

Using symmetry along ^ the integration on the periodic coordinate is performed directly 
and we find 

^■^This relation can be considered as some kind of generalized Smarr-like relation. As it was discussed 
in [12j . in the case of nonlinear electrodynamics (at least simple) Smarr relations do not exist. 
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/buik = -^/ i^dW-^l A\^){2XdxL-L)i^dW. (75) 
At this stage we should make use of the relation (I73|) which gives 

b^ = }-MB-TS- (76) 

up to terms of order 0{r^^). Direct computations of the other Euclidean terms give 
the following results 



^ = 2M-r^--^, (77) 



h^ = -M + rs + "f, (78) 

^ = ^nQ, (79) 

up to terms of order of 0(r^^). So for the Euclidean action, up to terms of order of 
0(r^^), we find 

4 = -MB + -M-i^-T5. (80) 
/3 2 2 2 £2 ^ ^ 

Finally, making use of f l65p we obtain 

^- = M-TS + 0{—) (81) 

which shows that in the limit — )■ oo we indeed have J//3 = M — TS. Let us note 
that the above considerations hold for arbitrary nonlinear electrodynamics not only for 
the Born-Infeld electrodynamics. 



C Off-shell considerations 

Let us consider different black holes, with different masses and entropie^. placed in 
a thermostat where the role of the thermostat is played by the Hawking radiation and 
the appropriate boundary conditions. Then we could ask which black holes would be 

-"^"^In the field black-hole thermodynamics it is usually accepted to work with the radius of the event 
horizon instead of the entropy. 
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in equilibrium with the thermostat. The equihbrium black holes, or in other words the 
equilibrium states are those for which the off-shel0 free energy 

Foff = M - T,^S, (82) 

has an extremum. In the above expression Tth is the temperature of the thermostat, it 
is kept fixed, and S and M are, respectively, the entropy and the mass of a black hole 
with radius of the event horizon vh- The off- shell free energy for several values of Tth 
is given in Figs. [21] and l22l 




Figure 21: The free energy F for sequence of black hole solution and the off-shell free 
energy F^s for several different temperatures of the thermostat Tfh as functions of the 
radius of the event horizon th- 

The equilibrium black holes have temperature of the horizon T that is equal to the 
temperature of the thermostat Tth. The different extrema are, actually, the different 
phases of the system. The free energy of the different phases, namely, the equilibrium 
black holes, is plotted with a thick line in Figs. [21] and [22] and it intersects Toff in its 
extrema. The stable phases are in the minima and the unstable in the maxima of Toff. 
The system is in the phase that realizes a global minimum of Tog. A phase transition 
occurs when with the varying of the temperature of the thermostat Tth the system 
passes from one minimum to another. The different branches of the equilibrium free 
energy T as a function of the temperature of the event horizon T (see, for example, 
the right panel of Fig. [T7|) can be obtained in the following way. As we know, the 
thermodynamic parameters, in our case the temperature and the radius of the black 
hole, are independent variables. They are related only at equilibrium, i.e. at an 

^^This energy is actually the off-equilibrium energy. The term off-shell is used to emphasize that it 
is not the natural free energy connected with the solutions. The on-shell energy, the one related with 
the solutions, is defined in the same way but Tth is substituted with the temperature of the event 
horizon T, which is different for black holes with different th- 
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Figure 22: Magnification of tlie enclosed region in Fig. |2T] 



extremum of the off-sliell free energy. The condition for an extremum of F^s is 



dr 



H 



= (83) 

rth=const 



This relation defines as an implicit function of Tth which, in the general case, may 

have several branches rH,a{Tth), = 0,1,2 When we take into account that at 

equilibrium T = Tth we obtain the th — T diagram (see, for example, the left panel of 
Fig. [IZl)- As we can see, the number of the extrema of Fofr for a given fixed value of 
T coincides with the number of branches of the th — T diagram for T = Tth and the 
positions of the extrema coincide with the radii of the event horizon of the different 
branches. 

The different branches of the on-shell or the equilibrium free energy (see, for exam- 
ple, the right panel of Fig. [T7|) can be obtained from the off-shell energy in the following 
way: 

Fa{T) = Toff (r^^,, (Tth), Ttt)\T,^=T ■ (§4) 

Let us consider, for example, black holes with P = 0.113. As can be seen in Fig. [21] 
for Tth = 0.173 the minimum corresponding to SBH and the one corresponding to 
VLBH are equal so this is the temperature of the phase transition Tph. 

The off-shell formalism can give also a very transparent description of the zeroth- 
order phase transition that occurs in the subinterval of the magnetic charge P^^ < P < 

,(2) 



Pp^ ■ The off-shell free energy is plotted in Fig. [231 The right panel is a magnification 
of the enclosed region in the left panel of the same figure. The three curves correspond 
to temperatures belonging, respectively, to the three intervals Tq < T < T^i^, T^in < 
T < Tph and Tph < T. It can be seen that for Tq < T < Triin the off-shell free energy has 
only two extrema - one minimum and one maximum [^. Hence, in that temperature 



16 



Notice that the on-shell free energy intersects the curve of the ofF-shell free energy only at two 
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interval only one locally stable phase is present - the VLBH. For temperatures in the 
interval Tmin < T < Tph a second minimum of the off-shell free energy corresponding to 
SBH occurs. It is lower than the minimum corresponding to the VLBH so the SBH is 
thermodynamically favored. As we have already mentioned the phase transition from 
the VLBH to SBH that occurs at T = Tmin is of zeroth order. As can be seen in Fig. 
[23| in the last temperature interval Tph < T the ensemble is again dominated by the 
VLBH since its minimum is lower than the minimum corresponding to the SBH and 
at the point T = Tph first order phase transition occurs. 




r 

;; 



Figure 23: The off-shell free energy for the case of the zeroth-order phase transition. 
The right panel is a magnification of the enclosed region of the left panel of the figure. 
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